Neutrino oscillations in a Robertson- Walker Universe with space time foam 
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In [lj we have studied decoherence models for flavour oscillations in four-dimensional stochastically 
fluctuating space times and discussed briefly the sensitivity of current terrestrial and astrophysical 
neutrino experiments to such models. In this addendum we extend these results to incorporate 
the effects due to the expansion of the Universe, so that our analysis can be useful in studies of 
extragalactic high-energy neutrinos, such as those coming from Gamma Ray Bursts at cosmological 
distances. Unfortunately for some microscopic models of foam, constructed in the string theory 
framework, we arrive at pessimistic conclusions about the detectability of the decoherence effects 
via flavour oscillation meaurements. 
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In [l[ we have discussed the propagation of both 
scalar and Dirac (or Majorana) flavoured particles in 
stochastically-fluctuating space times ("foam") 0], and 
we derived the corresponding expressions for the flavour- 
oscillation probabilities. The presence of a quantum foam 
results in modifications to the oscillation-probability for- 
mula, through damping factors in front of the oscillatory 
terms and alterations of the oscillation period. Their 
precise form (as functions of the energy of the probe) 
and strength depend crucially on the details of the mi- 
croscopic model used to describe the space-time foamy 
backgrounds. The non-observation so far of such effects 
in various neutrino experiments, both terrestrial and as- 
trophysical, places stringent constraints on the foam pa- 
rameters, which in turn leads to the exclusion already of 
some models. The majority, however, of the oscillation 
models of space time foam, involve damping exponents 
which are proportional to the square of the square-mass 
differences between the various flavours ((Am 2 ) 2 ), and 
as such is very suppressed. 

Sensitivity to such small effects could only come from 
cosmic neutrinos, which travel huge distances, of cosmo- 
logical scale, from their emission at remote celestial ob- 
jects till observation. Examples of such cosmic probes 
of quantum-foam decoherence can be provided by high- 
energy neutrinos from Gamma-Ray-Bursters (GRB) [3] 
at high redshifts, z > 1. Currently there are many phe- 
nomenological studies of such decoherence effects in up- 
coming high-energy cosmic neutrino facilities 0] , however 
the models used are simplistic and quite generic models of 
Lindblad decoherence [5| , without any attempt to discuss 
microscopic situations. Such an attempt was made in [l|], 
where the various damping exponents have been derived 
for some stochastic space-time backgrounds coming from 



microscopic models of space-time foam, including string- 
inspired ones 0. However, there the Universe expansion 
has been ignored. 

In the cases of cosmic probes the effects of the ex- 
pansion of the Universe cannot be neglected, and hence 
the analysis of [lj needs to be extended in order to in- 
corporate them. This is done in this brief note. We 
shall discuss first the case of scalar probes, deriving 
the decoherence-induced modifications to the oscillation 
probability in a Friedman-Robertson- Walker Universe, 
with space-time foam. As we shall demonstrate, if one 
considers the case of a slow-expansion rate for the Uni- 
verse, compared to the momentum/energy of the probe, 
which is a satisfactory approximation for high-energy 
probes in late eras of the Universe, such as the ones of 
interest to phenomenology, then the results for the eval- 
uation of the decoherence exponents for scalar fields can 
be carried through identical to the Dirac (or Majorana) 
fermion case. 



We commence our discussion by reviewing the differ- 
ent steps leading to the solution of the equation of motion 
for a scalar field in a spatially flat Robertson- Walker Uni- 
verse [fj , which is dictated by the astrophysical data and 
considered here. Expressed in terms of the cosmic time, 
the metric is 



g^ = diag(l, -a 2 (t),-a 2 (t),-a 2 (t)), 



The conformal time r/ is defined as 



g»v = c(v)v^, c(v) = a 2 (t), dn = 

a(t) 



(1) 



(2) 



The curvature scalar is R 



3c 3_ (6) 

2c V c) 
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The equation of motion for a scalar field <j> in this back- 
ground reads, 



V9 



9^ iVaa^duc, 



+ m 2 = -c 

c 



+ 4> - V 4> + cm 2 <t) = 0, (3) 



where g = | det («7 Mt -) | and the prime represents derivative 
with respect to the conformal time rj. We consider a field 
depending only on rj and x = x 1 and take the following 
ansatz 



<t>{ri,x) 



t>oe 



/c(?7) 



■Xk(v), 



(4) 



where 0o is a constant and k is the momentum of the par- 
ticle. With this normalization, no "friction term" Xk is 
generated from eq. , and we obtain the following equa- 
tion of motion for Xk 

Xk(v) + uUv)Xk(v) = 0, (5) 
where the time-dependent frequency is given by 



i 2 i 2 

k + cm 
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-if. 



(6) 



Note that, in the adiabatic approximation, where 
|c| 2 << fc 2 and |c| << /c 2 , particle production by the 
expanding Universe can be neglected. This approxima- 
tion, equivalent to neglecting the curvature R, is phys- 
ically reasonable for high-energy probes at late eras of 
the Universe, which is of phcnomenological interest. We 
then look for a solution of eq.© of the form 



Xk(v) = Fk(v) exp ( i / dr{ W k (j)') 



From eq.©, we then obtain: 
F + 2iFW + F{iW 



W 2 +LU 2 ) = 0. 



(7) 



(8) 



The latter equation is solved in the framework of the 
WKB approximation, for which the small parameter is 



(9) 



In this expansion framework, the various orders are de- 
fined by the order of conformal-time derivatives of F and 
W. To lowest order in the WKB approximation, we ne- 
glect F in eq.©, which leads to: 



• terms without derivatives: 
that Wk(rj) = Wfc(ry); 



-W 2 + u? 



0, such 



• terms with one derivative: 2FW + FW = 0, such 
that F k {rj) = l/y/w k (fi). 

The solution for the initial scalar field to lowest order in 
the WKB approximation, reads: 



<t>{ri,x) 



t>oe 



V c (v)^k(v) 



exp 



i / dr]' uJk(v') 



(10) 



where t]q denotes the emission time of the probe. We 
express the solution (fl"0|) in terms of the redshift z, which 
is related to the scale factor a(t) by the usual relation 



a(t) = = 



1 + z' 



(11) 



where a(now) = 1, when z = 0. In terms of the Hubble 
parameter H(z) = a/a, we obtain to lowest order in the 
WKB approximation, 

= (1 + z)+f^ exp f-i f Z dz'^1) , (12) 

where z$ > z (for 770 < T]), with the suffix here denot- 
ing quantities at emission, not to be confused with the 
(opposite) notation used in astrophysics. In our a nota- 
tion, the normalization is that at the time of observation 

2 = 0. 

To study the effect of the stochastically fluctuating ex- 
panding Universe on flavour oscillations, we consider for 
brevity the case of two flavours \<fr a > (a = 1, 2), related 
to the energy eigenstates \ fi > (i = 1, 2) by 

\<p a >= u «Ah > < <t><*\ = J2 U ^ < M 

i i 

\fi >=J2 U ^ > < fi\ = E^« < m -( 13 ) 

The unitary matrix U corresponds to the tranformation 
from the flavour eigenstates basis to the energy eigen- 
states basis. The states fi > have masses m, and satisfy 
the equation of motion © . The corresponding solutions 
for the frequency and the amplitude are denoted Wi{rf) 
and Fi{rf). If we assume the following pure state density 
matrix for the initial condition (no summation over a): 



Pivo) = \<t>a(m) >< <f>a(vo)\ 

= J^t/criE^I/iM >< /i(r?o) 



(14) 



the time evolution of the energy eigenstates leads then to 
the following evolution for the normalized density matrix 

pin) = X>ait/^|Mr?) ><fM\ 

ij 

- STtt tt* If (r, \ ^ f C« ^l c ^°) Fi (n)Fj(ri) 



c( V ) Fi^Fjirjo) 



xexp (i / dn'iWiW-Wjir)')) 



E u ai u* aj u* 0i u^\Mvo) >< Mvo)\ ^ fif^pK 

f^ 3 <v) Fi{Vo)FAno) 



xexp J drj'iWi^-Wjiv')) 



(15) 



As a consequence, starting from the initial flavour a, the 
transition probability to the flavour (3 ^ a is 



"(.Vo) 
c(rj) 



FMFM 
Fi(vo)Fj(r)o) 



xexp i dr)'(Wi(r)') - Wj(r,')) (16) 



Note that, if rj = 170, we obtain the standard result of 
flavour oscillations, 



Pn^P = UaiUajUpiUfij = 0, 



(17) 
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which can be readily verified, for example, with the choice 
/cosS -sin (A 

V sinS cos 6* J v ' 

Also, for any rj, equality of masses m, = rrij leads to 

Pa -" = cwmm) 5 UaiU ^ u ^ = °- (19) 

We consider now the effect of a stochastic background 
by introducing the following modification of the inverse 
metric, in terms of the conformal time, 



tin — 1/ \~uv ~ uiy niy . i iiv 



(20) 



and the stochastic components are coordinate- 

independent and are assumed to be small compared to 
unity. The equation of motion is then 



gOO [j> + -J>j + g 01 \ 2d x <t> + ^d x <f>) + g 11 d 2 <f> + c( V )m 2 <f) = 0. 

(21) 

The ansatz (Q} leads to 



g 00 X k(v) + 2ikg (n Xk (v) + tlUriXkiri) = 0, (22) 



where 



mv) = -gk 2 + c(r,)m A ~g m ^R 

D 

2/ i ,11 , 2 .00 c r, 

= uj k (rj) -h k -h -R. 



(23) 

We look for a solution of eq. (|22p of the form 

Xkirf) = F k (r,)exp (i £ drf W k ( V '^J , (24) 
such that F and W satisfy 

-oo + 2i p w + iFl y _ FW 2^ +2ikg -°i (p + lFW ^ + tf F = 

(25) 

In the framework of the WKB approximation, neglecting 
F leads to: 

• terms without derivatives: — g 00 W 2 — 2kg 01 W + 
Q 2 = 0, such that 

%01 i 



W = " fc foo + ~oo \A 2 (s 01 ) 2 +ff 00 « 2 ; (26) 



• terms with one derivative: g 00 (2FW + FW) 
2kg m F = 0, such that 

1 



F = 



yjkg 01 + g 00 W 



(27) 



If we consider only the linear terms in the fluctuations 
W v , the latter expressions reduce to 



2 2a, 2 fa) 



kh u 



„, . 1 Aft 00 h^k 2 



(28) 



The time-dependent part of the probability (Tl6[) is 

(29) 

and should be integrated over the random fluctuations of 
h^ v . For concreteness, in this note we consider a Gaus- 
sian distribution for the latter. Since only h 00 and h 11 
appear in the expression (|29p . we define the average 

<•••>= J- / dh 00 dh}\---)ex V (-{h w /af-{h lx /af). 
irc^ J 

(30) 

Note that, in order to define the distribution, we con- 
sider quadratic terms in W whereas we kept only the 
first order in in the result (f2"5|) , but this is not con- 
tradictory if a is of order h^ v . The lowest order of the 
WKB approximation gives then 



< P v > 



c(:no) 

cfa) 



exp^ij dr)'(ui(r)') -W;,fa'))j 



X 1- i- 



.k 4 a 2 



Ayfa) [ drj' (l/u;i(ri') - lMfa') 



x exp 



16 



drj' (uii(r)') -Wj(rj') 



a 2 k' 
16 



^"^'(im^O-iM^'^'msi) 



where 



Ay fa) = 



1 1 

+ ■ 



<<fao) <fa) u;?fa ) a; 2 fa) 



(32) 



In order to have a more intuitive understanding of this 
result, we neglect the curvature R and take cu 2 ~ k 2 + 
cm 2 . We make then an expansion in |m, — rrij\/k and 
obtain, in terms of the redshift, 



< P z >= 



1 + 2 



2 2 
m i - m i r, \ 

exp <; j — — 9 J(z ; 



1 + 2o y " _jr 1 " 2/t 2 

2 __4 4 2 _2 



(33) 



j j i * m i ~ m - j, z s 4 - 2 + 2(2o - 2) 



16 fc 4 



(l + Zo ) 2 (l + 2) 2 



a 2 / m 2 - m 2 \ 2 r2 



l + z \ 2 / a 2 /m 2 _m|\ 2 2 



nij — rrij a rrq — rrii z 2 — z + 2(z — z) 



x exp 1 
where 



2fe 2 



16 fc 4 (l + z ) 2 (l + z) 2 



7(2) 



7 W = fc y^*/c(,) = -* yao H(y)(1 + ^ )a 



(34) 



zq (z) denotes the redshift at emission (observation) of 
the particle and we have assumed that a 2 <C 1, which, as 
we shall see later, characterises semi-realistic microscopic 
models of foam. The oscillation probability is, of course, 
the real part of the above expression, and thus we ob- 
serve that the effects of the foam, apart from introducing 



4 



damping factors, also modify the oscillation period (c.f. 
imaginary exponents). Notice also that a 2 is allowed in 
general to depend on the energy of the probe, which is 
found explicitly in some string- inspired models @, B • 
In ([33| , the first line corresponds to the oscillation term 
in a smooth expanding Universe; the second line leads 
to a correction to the phase of these oscillations, arising 
from the metric fluctuations, and the third line represents 
a damping, also generated by the random fluctuations of 
the metric. Note that, for a flat space time, z = zq 
and I(z) = k(t — to), such that we recover the results of 
[l( . Assuming the standard Cosmological-Constant-Dark 
Matter Model (ACDM), 



H(z) = # \/fiA + (l + ^ 3 ^ 



(35) 



where the various energy density parameters are de- 
fined, as usual, in units of the critical density of the Uni- 
verse, we thus obtain for the case of flavoured particles, 
such as neutrinos, emitted at a period in the Universe 
with red-shift zq, 



I(z) 



k 



da 



1 



(1 + z') 2 + (1 + z'^Qm 



(36) 



To extend the above result to the case of neutrinos 
(Dirac or Majorana) we first note that the motion of 
fermions in a curved space time is described in terms of 
the vierbeins e^, which span the flat tangent space time 
at each point x^ 1 of the manifold, and are related to the 
metric by g^ v = e®e^r) a b, where Latin indices refer to 
the local inertial (flat, tangent) frame and are contracted 
with the Minkowski metric r\ a b, whereas Greek indices are 
contracted with the metric g^ v . The equation of motion 
for fermions in a curved space time is [6| 



i7 a e£ (<9 M + T M ) ip - mip = 0, 
where the spin connection is 



(37) 



(38) 



and the gamma matrices are defined in the local inertial 
frame, and therefore satisfy 



{7 a ,7 6 } = 2r? a 



(39) 



With the conformally flat metric ([20)1 . it is easy to check 
that 



+ 0{h 2 ). 



(40) 



If we neglect the time derivatives of the scale factor, the 
Christoffel symbols vanish, as well as the spin connec- 
tions. The equation of motion for fermions, then, reads 
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d a + -Kd, 



o. 



(41) 



and a multiplication by the complex conjugate operator 
leads to 

(d"d P + hTdpdv + c(r;)m 2 )v = {sTd^du + c(r?)m 2 )v = 0. 

(42) 



This equation is similar to eq. plj) for a scalar field, when 
the derivatives of the scale factor are neglected when 
compared with the momenta/energies of the probe, and 
hence in this approximation of slow expansion, the pre- 
vious results for the decoherence damping exponents of 
the scalar field carry through intact to the fermion case 
(above we have discussed for concreteness Dirac fermions, 
but the extension to the Majorana case is straightforward 
and does not alter the final result). 

Let us now make a few remarks on the result (|33p . 
([36]) . As charactertistic of the Gaussian space-time fluc- 
tuations [l|, |2j, the damping exponent depends on the 
square of the time traversed by the probe, i.e. the damp- 
ing terms are of the form exp \—D t 2 ). Should one use 
other distributions for the space-time fluctuations, for 
instance the Cauchy-Lorentz distribution that could 
arise in some string inspired models of quantum foam [2| , 
then the decoherence damping assumes a Lindblad ex- 
ponential form [5j. From (|33[) . the decoherence-induced 
damping exponent is of the form: 



exponent = 

1 32 



A:-' 



k 



dz 



H(z') (1 + z') 2 



(43) 

Some remarks are in order concerning the nature of the 
variance a 2 . Its precise estimate requires microscopic 
models of the stochastically fluctuating space time. 

One such model is inspired from string theory, namely 
the D-particle foam model [7|, in which a space time 
fluctuation arises from topologically non-trivial interac- 
tions of open strings (representing matter excitations on 
a brane world scenario) with D-particle defects. The sce- 
nario involves capture and re-emission of the open string 
by the defect, which is microscopically represented by 
the appearance of intermediate string states stretched be- 
tween the D-particle and the brane world. This leads [8( 
to causality-respectful time delays, proportional to the 
energy E of the incident string, 



At ~ dE ■■ 



E 



(44) 



On the other hand, the recoil of the D-particle defect 
during the above process implies local distortions of the 
space time, which are of Finsler type, in the sense of 
depending on both the coordinates and the momentum 

transfer Afc of the string state. In particular, one obtains 
non diagonal metrics, which for small recoil velocities cor- 
respond to deviations from the Minkowski metric of the 
form 



% = 1,2,3 



(45) 



As discussed in detail in [9(, quantum fluctuations of 
the recoil velocity arise by summing up over world-sheet 
genera, which translate in turn to stochastic quantum 
fluctuations in the induced finsler space time. The cor- 
responding distributions are Gaussian, with a variance 
a 2 , which, over the time scale of the interaction of the 
D-particle with the string state St (T4"4")) . becomes 



cry = (UiUj) = a d i:j = g 



.(At)' 



c _ 2 E x 

~ 9s M 2 



(46) 
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Using (|46j) in (|43|1 . then, and approximating to lead- 
ing order in small quantities E ~ |fc (which for the case 
of high-energy neutrinos is a perfectly good approxima- 
tion), we obtain for the decoherence damping exponent 
the estimate: 



exponent 



32MSH, 
dz 



m i) X 



y/Q A + (1 + z') s Q M (1 + ^') 2 



(47) 



for the ACDM cosmology, where we set the red-shift of 
observations 2 = 0. 

We first observe that for this model of space-time foam, 
the decoherence damping exponent is independent of the 
probe's energy or momenta. Taking into account that 
the value of Hubble parameter today Hq ~ 10~ 24 Mp, 
where Mp ~ 10 19 GeV is the four dimensional Planck 



mass scale, we can express (|47| as: 



io 48 (ff 



(Am 2 ) 2 
Mi 



For g s Mp/M s ~ 1, which is a phcnomcnologically real- 
istic situation, allowing also for large string scales and 
very weak string couplings, and taking into account that 
the fc-independent red-shift factor in the last line of Wf\ 
is of order one for red-shifts zq < 5 we are interested in 
here, one observes that the resulting decoherence coeffi- 
cient is unobservably small (< 10 -72 ), for phenomenolog- 
ically realistic values of mass-squared differences between 
neutrino flavours of order Am 2 ~ 10~ 3 — 1CU 5 eV 2 . In 
fact to obtain decoherence exponents of order one, one 
needs probes that have travelled the age of the observ- 
able Universe. This has to do with the fact that the factor 
^ A ^d is numerically of the same order as the (observed) 

Cosmological constant 1CU 122 MJ, ! In a similar way the 
other decoherence corrections in (|3"3")) are also unobserv- 
ably small for this model. Thus, this model of D-particle 
foam cannot be falsified via decoherence damping tests. 

On the other hand, for other distributions of recoil 
velocities, such as the Cauchy-Lorentz (CL) distribution, 
with characteristic function <&(£) = -^^ppy, used in 

one has exponents of the form: 



Am 2 

exponents ~ -yH^ 1 -^- (I(z)H /k) (48) 

that is, exponents of Lindblad type, varying linearly with 
the (cosmic) time. Again the fc-independent red-shift 
time factor I(z)Ho/k (c.f. (j3"o| ) is contributing terms 
of order one for the region of red-shifts of interest z < 5. 
The CL type distributions may be encountered [l[ in pop- 
ulations of D-particle defects, as opposed to the virtual 
Gaussian quantum fluctuations we have discussed previ- 
ously. 



Since the parameter 7 is essentially arbitrary in such 
models, one might hope that the CL distributions may 
lead to observable effects for cosmological distances 1 < 
z < 5, for infrared neutrinos. In fact, for 5m 2 £ {10~ 3 — 
1CT 5 } eV 2 = 1CT 59 - 1CT 61 M 2 , H ~ 10" 24 M P we ob- 
tain the following estimate for the decoherence damping 
exponent of the CL distributions, modulo the order-one 
I{z)H„/k factors: ^(1(T 35 - 1CT 37 ) . This is of order 
one for neutrinos of eV energies, if 7 <~ 10 7 for the lower 
bound, which however is an unnatural value for a disper- 
sion on velocities. For values 7 < 1, becomes of order 
10~ 4 for the lower bound, only for infrared neutrinos of 
energies in the meV or less. 
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